Amorphous Precursors of Crystallization during Spinodal Decomposition 
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A general Landau's free energy functional is used to study the dynamics of crystallization during 
liquid-solid Spinodal Decomposition (SD). The strong length scale selectivity imposed during the 
early stage of SD induces the appearance of small precursors for crystallization with icosahedral 
order. These precursors grow in densely packed clusters of tetrahedra through the addition of new 
particles. As the average size of the amorphous nuclei becomes large enough to reduce geometric 
frustration, crystalline particles with a body center cubic symmetry (bcc) heterogeneously nucleates 
on the growing clusters. The volume fraction of the crystalline phase is strongly dependant on the 
depth of quench. At deep quenches, the SD mechanism produces amorphous structures arranged in 
dense polytetrahedral aggregates. 



According to the classical theory of symmetry breaking 
phase transitions, the crystallization process may pro- 
ceed via two distinct mechanisms, nucleation and growth 
(NG) or spinodal decomposition (SD) In the NG 
mechanism the initial disordered phase is metastable and 
the relaxation of the system is promoted by the overcom- 
ing of a free energy barrier. This process involves the 
formation of a critical nucleus of the crystalline phase 
by structural fluctuations. Differently from the NG pro- 
cess, SD does not require large fluctuations to initiate 
the phase transition and it is characterized by the ex- 
ponential growth of density fluctuations of a dominating 
wavelength, entirely determined by the thermodynamic 
properties of the system. At early times this process 
leads to the formation of disordered isotropic states that 
eventually evolve towards the equilibrium phase through 
coarsening HI- 

Although this classical picture provides a qualitative 
description of the phenomenology involved in the phase 
transition process, during the last thirty years, experi- 
ments, theory and simulations have shown that the pro- 
cess of crystallization is far from being trivial. For ex- 
ample, in the NG regime, the kinetic pathway towards 
equilibrium may involve intermediate phases with a sym- 
metry different than the corresponding to equilibrium [3]- 
P|. According to the Ostwald's step rule the nucleated 
phase is not necessarily the thermodynamically most sta- 
ble, but the energetically closest to the disordered state 
0. A similar process was observed in the crystalliza- 
tion of globular proteins and colloids, where density in- 
homogeneities in the fluid lower the free energy barrier 
for NG On the other hand, the relaxation in the 

spinodal region is still poorly understood. For example, 
it has been recently observed that thermal fluctuations 
play an important role in the dynamics of SD. In the 
neighborhood of the spinodal line it was found that the 
equilibrium phase can be pseudo-nucleated by density 



wave fluctuations [10]. A complex phase separation ki- 
netics have been observed also in different systems, like 
colloids and polymers, where it has been found that a 
spinodal-like dynamics can precedes the crystallization 
[if 
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process 

In this work a continuous Landau's free energy expan- 
sion and a relaxational dynamics are used to study the 
kinetics of crystallization of a bcc structure in the spin- 
odal region. As compared with molecular dynamics sim- 
ulations, the phase field approach employed here natu- 
rally incorporates the elasticity of the bcc crystals and 
provides an efficient approach over diffusive time scales, 
while the identification of the precursors of crystalliza- 
tion as well as amorphous and crystalline regions can 
be easily determined through the local amplitude of the 
continuum order parameter. 

Symmetry-breaking phase transition in a wide variety 
of systems can be studied through an expansion of the 
free energy $ in terms of an appropriated order parame- 
ter tp. Here we employ the expansion [12j : 



^(r)V(r') 



dV'}dV, (1) 



where W{ip) = —ri/j'^ + i/ip^ + Xip'^. Here the order param- 
eter represents the deviation of the particle density from 
the uniform value characteristic of the liquid state (where 
-0 = 0). The parameter r is proportional to T5 — T, be- 
ing a measurement of the deep of quench, with Ts the 
spinodal temperature where the continuous phase tran- 
sition begins. The constants v and A are related to the 
symmetry and saturation of ip at equilibrium, and 13 is a 
free energy penalization to spatial variations of the order 
parameter 0) 01, [isj - 

The dynamics of SD can be studied through a Cahn- 
HiUiard equation [T^: dtji/dt = AfV2{(5$/5V}, where M 
is a mobility coefficient [l5|. 
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Since the initial disordered state is characterized by 
small fluctuations {ip ~ 0), the early dynamics is al- 
most linear and the system's state can be described as a 
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FIG. 1. Early autocorrelation function C(r) for a system 
quenched into the spinodal region with reduced temperatures 
Tr = (t - ts)/ts = 1x10"'^ (black line) and = 1x10"^ 
(blue line). The dashed red line indicates C(r) ~ 1/r. The 
inset shows the early time real space fluctuations (left) and 
the formation of precursor by the non-linear enhancement of 
density fluctuations (right). 



random superposition of density waves of the form [lOj: 
V'Cr, t) = A^exp {ik - r + X{k) t). Here is the ini- 
tial amplitude of the k-mode and the ampHfication fac- 
tor X{k) — —D k* + Tk^ — b selects the range of unstable 
modes (those modes for which A(fc) > 0). In real space 
the system display a disordered pattern characterized by 
dominating length scale related with the most unstable 
modes p^ . 



For systems quenched near below the spinodal line, 
where ts = 2\/bD, the region of unstable modes be- 
comes sharply peaked around a dominating wave vector 
amplitude ko ^ \J t / D. In this case the random super- 
position of modes leads to the emergence of a strongly 
correlated filamentary network of density wave fluctua- 
tions, with strong similarities to those found in quantum 
billiards and other physical phenomena involving random 
wave superposition [lO[. Figure 1 shows the long range 
density correlations that emerge in the early state. In 
this case, the azimuthally averaged two point correla- 
tion function, C{r) = (J dr'-0(r')-0(r+r')), of a critically 
quenched system behaves like a Bessel function decaying 
as C(r) 1/r, in agreement with the theoretical pre- 
dictions for the random superposition of waves • For 
sub-critical systems (r > Tj) the network of fluctuations 
loose correlation and C(r) decays faster. As shown be- 
low, the early correlations in ij: have a profound effect in 
the later evolution of the system. 

As time proceeds, there is a continuous amplification 
of il) until the anharmonic terms of the free energy func- 
tional cannot be neglected and nonlinear dynamics comes 



into play. While at deep quenches there is a lack of corre- 
lation in ip and the classical picture of SD is recovered, at 
shallow quenches the early network of density wave fluc- 
tuations trigger the inhomogeneous appearance of pre- 
cursors for crystallization (Figure 1). Note that the local 
symmetry of the precursors is dictated by non-linear dy- 
namic effects and does not necessarily coincides with the 
symmetry of the phase of equilibrium. 

In order to track the evolution of the system we iden- 
tify structural features by means of Voronoi diagrams 
through the centers of the particles (bulk particles are 
defined as local maxima of Voronoi diagrams indi- 




FIG. 2. (a) Precursor formed by four icosahedral particles 
(red) and their first neighbors (yellow). Red lines indicate 
bonds between icosahedral particles and blue lines between 
icosahedra and first neighbors. Inset: Voronoi diagram of an 
icosahedral particle (dodecahedron), (b) Precursor formed 
by a compact cluster of icosahedra arranged in the vertex 
of a tetrahedron, (c) Propagating amorphous nuclei (yellow) 
having compact clusters and small chains of icosahedral par- 
ticles (red), (d) Regular tetrahedron (down right), where yel- 
low spheres indicate the particles and green sphere indicates 
the tetrahedron center's, and polytetrahedral aggregate (up), 
where black lines connect tetrahedra with common faces. 
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cate that the early precursors are constituted by com- 
pact clusters and chains of icosahedrally arranged par- 
ticles (Figs. 2a y 2b). That is, the symmetry of the 
precursors is dictated by the non-linear dynamics rather 
than by equilibrium states. 

The presence of precursors with icosahedral symmetry 
in the early stage of SD is not surprising. One of the old- 
est scenarios for the formation of amorphous matter and 
glasses is geometrical frustration. According to this sce- 
nario, the amorphization of supercooled liquids is related 
with their tendency to form locally compact structures, 
like icosahedra, incompatible with translational symme- 
tries The symmetry found in the precursors is 
also in agreement with the Alexander-McTague predic- 
tion which indicates that local icosahedral order is fa- 
vored in the absence of a NG mechanism . 

Once the precursors have been formed, they begin to 
grow with an approximately spherical shape at a constant 
rate by the addition of new particles (Fig. 2c). Given 
that the initial symmetry is incompatible with transla- 
tional order, the aggregation process produces growing 
amorphous clusters without any appreciable symmetry. 
As the size of the precursors increases, the geometrical 
frustration is reduced and a number of particles begin to 
aggregate with a crystalline structure at the front of the 
propagating precursor (Fig. 3a) j^lj . 



This process of surface crystalline nucleation shows 
similarities with Monte Carlo results for the heteroge- 
neous nucleation on the surface of spherical colloidal 
seeds [22|. In agrement with the Monte Carlo simula- 
tions, here we also observe that the crystalline particles 
do not span the whole surface of the precursor because 
crystals cannot grow without generating amorphous re- 
gions. Then, the growing nuclei shows an amorphous 
core surrounded by a shell constituted by amorphous 
and crystalline phases (Fig. 3a). A similar phenomena 
has been also observed in the freezing gold nanoclusters 
where initially the system crystallize on the surface of 
the nuclei while the amorphous core with icosahedrally 
arranged particles remains stable for some time 
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Figure 3b shows the average radial profile of the ratio 
between the fraction of crystalline particles on the prop- 
agating nuclei /c^ and the final fraction of crystalline 
particles throughout the system fcrp. This figure shows 
that crystal particles begin to nucleate at the surface 
of the amorphous precursors for radius of the order of 
'^c ~ (3 — 5)a, with a the inter-particle average distance, 
also in good agreement with the results of Cacciuto et al. 
[23 [ . Then for r < Tc the high curvature of the precursors 
clearly frustrates the formation of crystal bonds. The 
subsequent growth and collision of the different propa- 
gating nuclei leads to the formation of a structure having 
crystalline and amorphous regions (Fig. 3c). 

Once the SD process has been completed, the volume 
fraction of crystals in the system depends strongly on 
temperature. As the depth of quench increases, there is a 




FIG. 3. (a) Propagating nucleus formed by an amorphous 
core (yellow particles) and crystalline particles (blue parti- 
cles), (b) Typical nuclei density profile, (c) System com- 
posed by amorphous and crystalline regions after the growth 
and collision of the different nuclei, (d) Fraction of crystalline 
particles as a function of Tr- (f)-(g) <16 and distributions 
for the amorphous (continuous line) and crystalline (dotted 
line) regions. 



larger number of unstable modes in the system and C(r) 
decays faster than 1/r (Fig. 1). Consequently, as the 
temperature of quench drops, the number of amorphous 
precursors increases, incrementing in this way the frac- 
tion of amorphous material. Figure 3d shows the fraction 
of crystalline particles /cy as a function of the reduced 
temperature — (t — Ts)/ts. We found that fcrp fol- 
lows a power law with (/cy ~ ). This behavior 

is consistent with previous results indicating that the av- 
erage distance between precursors ^ diverge as ^ ^ 
{r, ^ 1/5, fcr^e) \M- 

To identify the symmetry of the crystalline phase sur- 
rounding the amorphous precursor we applied the bond 
order parameter analysis introduced by Steinhardt et al 
where the local structure is characterized in terms 
of the symmetry of near neighbor bonds by using spher- 
ical harmonics. The distributions of the bond order pa- 
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FIG. 4. Pair correlation functions g{r) calculated with the 
centres of all the particles (red) and with the centres of the 
icosahedra (black). Inset: distribution for the qs bond order 
parameter (left) and the pair correlation function calculated 
through the center of the regular tetrahedra (right). 

rameters are very sensitive to the underlaying symmetry, 
allowing a clear identification of the structure. Figures 
3f and 3g show the distributions of the rotational invari- 
ants qe and wq , for the amorphous and crystalline re- 
gions. These distributions reveal that crystal particles ar- 
range in a bcc lattice (peaks at ~ 0.5 and wq ^ 10~^) 
while in the amorphous phase both distributions become 
broader and shifted towards the low bond order param- 
eter region [1], [l7{ . 

At long times the relaxation of the system continues 
through a complex process involving structural changes 
(coarsening) , where there is an invasion of the crystalline 
structure on the amorphous regions. However, at low 
temperatures the dynamics can be highly arrested and 
the structure freezes. 



In order to study the glasslike amorphous structure in 
the deep spinodal region, here we also analyze the struc- 
tural features of a system where the formation of crystals 
is totally frustrated (r^ > 10~^ in Fig. 3). Figure 4 shows 
the pair correlation function g{r) for a typical amorphous 
structure. In this figure we also include the pair corre- 
lation function associated to those particles with icosa- 
hedral symmetry. Note that the amplitude of the first 
maximum of the pair correlation function calculated with 
the icosahedra is about three times larger than the am- 
plitude of the first maximum of g{r) calculated with all 
particles. In addition, the position of the first maximum 
for the icosahedral pair correlation function is slightly 
smaller than the corresponding to all particles, indicat- 
ing that the icosahedral bonds are under compression. 

The structural features found here (small clusters and 
chains of icosahedra and the behavior of pair correla- 



tion functions), agree with experimental findings in su- 
percooled colloidal systems undergoing NG ^] , jM] , sug- 
gesting that icosahedra could also play a role in the amor- 
phization of those systems. 

The structural characterization of amorphous matter 
has been one of the major problems in condensed matter 
physics during the last fifty years. In this sense, one of 
the best geometrical descriptions is based in aggregates of 
polytetrahedral clusters which have been invoked in order 
to describe the structure of liquids, glasses and quasicrys- 
tals. 

The geometrical features of the tetrahedra can be ob- 
tained through Voronoi diagrams (Fig. 2d). To test for 
the regularity of tetrahedral configurations in the amor- 
phous phase we calculate the q^ order parameter, which 
is sensitive of local tetrahedral order around the tetrahe- 
drons center (93 = \/5/3 for a regular tetrahedron). The 
inset of Fig. 4 shows the ga-distribution for our system. 
The sharp peak of the distribution at 53 ^ 0.73 clearly in- 
dicates that regular tetrahedra are the fundamental units 
of the amorphous states obtained via SD. The analysis 
of the tetrahedral configurations also shows that most of 
the particles (~ 90%) are involved in at least one regu- 
lar tetrahedron, indicating the presence of an amorphous 
|)lMise characterized by a dense polytetrahedral structure 

The inset of Fig. 4 also shows a pair correlation func- 
tion calculated through the centres of tetrahedra. After 
a few oscillations g{r) rapidly decays towards the asymp- 
totic value for r > a, showing that the tetrahedra aggre- 
gates in local configurations without long range order. 
By using real space plots we observe configurations of 
tetrahedra packed on common faces, forming disordered 
clusters (Fig. 2d). Similarly to dense disordered packing 
of hard spheres, here we also found that the aggregates 
are mainly formed by five-fold rings obtained by packing 
five tetrahedra around a common edge [25j . 

In conclusion, we have presented evidence of a new 
mechanism of crystallization during SD, induced by 
amorphous precursors. Local icosahedral clusters and 
amorphous precursors are kinetically favored by the early 
dynamics of SD. Crystalline particles aggregate on these 
precursors in a similar way as heterogeneous nucleation 
on spherical colloidal seeds. Our results are in qualita- 
tive agreement with recent NG work in colloidal systems 
[26| . Although the relaxational mechanism studied here 
can be easily confused with conventional NG, it is hoped 
that the phenomenon of crystallization in the neighbor- 
hood of the spinodal can be experimentally observed in 
systems with a slow dynamics and a high length-scale 
selectivity like hard-sphere colloidal suspension or block 
copolymers. 
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